UNCLASSIFIED 


AD  NUMBER 

AD232705 

NEW  LIMITATION  CHANGE 

TO 

Approved  for 
unlimited 

public  release,  distribution 

FROM 

No  Foreign 

AUTHORITY 

ONR  ltr.,  15  Jun  1977 


THIS  PAGE  IS  UNCLASSIFIED 


flruiBi)  Services  Technic*!  Infenealien  Agency 

ARMhGrON  HALL  STATION }  ARLINGTON  12  VIRGINIA 


NOTICE;  Wl! 
)raiRKATi 
*/rni  a  n 

THE  U.  L  A  T 
OBLIGA11C  N  V 
HAVE  FVRJ4I11 

~*Myvn.-8,  at 

IMPLICATE  y*f 
OR  ANV  r.  TK  V 
i’LXocU!  ■; 
THAT  M*i.  f  IK 


m  OOVIRNIONT  OR  OTHER  DRAWINGS,  BPECIFTCATiaNfl  OR 
A  IE  UBD  FOR  ANT  PURPOSE  OTHER  THAN  IN  CONNECTION 
[  RELATED  GOVERNMENT  PROCUREMENT  OPERATION, 

1 RNMENT  THEREBY  INCURS  NO  RESPONSIBILITY,  NOR  ANY 
1 1ATSOEVER;  AND  THE  FACT  THAT  THE  GOVERNMENT  MAY 
U  tTED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE  SAID 
f  CIFXCATEONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
Z  R  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER 
I  PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR 
C  MANUFACTURE,  USE  OR  SELL  ANY  PATENTED  INVENTION 
\  NY  WAY  BE  RELATED  THERETO. 


§3 

o 

o 


Interim 


Report 


REFINEMENTS  OF  THE  THEORY 
OF 

THE  INFINITELY-LONG,  SELF-ACTING, 
GAS- LUBRICATED  JOURNAL  BEARING 


I-A2049- 10 


c 


*  .  'Wi 


by 

Harold  G.  Elrod,  Jr. 
Albert  Burgdorfer 

January  1960 


l 

U 

-J 


* 

2 

< 

i 

V? 


*- 

«n 

< 


X 

O 

K- 

£ 


3 

£ 

o 

ft! 

> 

W 

Z 

o 

»- 

o 

* 

►J» 

ft* 

4./ 


W1 

v* 


Prepared  under 

Contract  Nonr-2342(00) 
Task  NR  097-343 


Jointly  Supported  hy 

DEPARTMENT  OF  DEFENSE 
ATOMIC  ENERGY  COMMISSION  J . 

MARITIME  ADMINISTRATION  !  ;  : 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 

i  1  *  ’ 

Administered  by  i.,'.  v  / 

OFFICE  OF  NAVAL  RESEARCH 
Department  of  the  Navy 

COPY  NO.  103 


THE  FRANKLIN  INSTITUTE 

LABORATORIES  FOR  RESEARCH  AND  DEVELOPMENT 
PHILADELPHIA  PENNSYLVANIA 


THE  FRANKLIN  INSTITUTE  *  Laboratories  for  Research  and  Development 


Interim  Report 


I-A2049-I0 


REFINEMENTS  OF  THE  THEORY 
OF 

THE  INFINITELY- LONG,  SELF-ACTING, 
GAS-LUBRICATED  JOURNAL  BEARING 

by 

Harold  G.  Elrod,  Jr. 

Albert  Burgdorfer 


January  i960 


Prepared  under 

Contract  Nonr-2342( 00) 
Task  NR  097-343 


Jointly  Supported  by 
DEPARTMENT  OF  DEFENSE 
ATOMIC  ENERGY  COMMISSION 
MARITIME  ADMINISTRATION 
NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 


Administered  by 

OFFICE  OF  NAVAL  RESEARCH 
Department  of  the  Navy 


Presented  at  1st  International  Symposium  on  Gas  Lubricated 
Bearings  at  Washington ,  D.C. 

October  26  to  October  28,  1959 


REPRODUCTION  IN  WOLE  OR  IN  PART  IS  PERMITTED  FOR  ANY  PURPOSE  OF  WE 

U.S.  GOVERNMENT 


THE  FRANKLIN  INSTITUTE  •  Laboratories  for  Research  and  Development 


I-A20i)9-10 


ABSTRACT 


The  lubrication  equations  for  an  arbitrary  Newtonian 
fluid  are  derived  directly  from  the  general  equations  for 
conservation  of  mass,  momentum,  and  energy.  From  the  lubri¬ 
cation  equations  an  inequality  is  obtained  for  the  internal 

film  temperature  rise.  The  isothermal  film  equations  are 
.  U'd 

then  derived.  Then,  for  perfectly-aligned  self-acting  jour- 
nal  bearings,  a  conservation  equation  is  obtained.  For  gas 
bearings  this  condition  gives:  . .  i 

i 

2  it 

f  P^h3d0  =  constant 

O 

along  the  axis  of  the  bearing.  Application  of  this  condition 
to  the  infinitely-long  gas  bearing  gives  more  accurate  pres¬ 
sure  solutions  for  this  case. 

The  Katto-Soda  form  of  the  differential  equation  for  the 
infinitely-long  bearing  is  solved  by  a  series  expansion  in 
the  eccentricity  ratio,  the  first  terms  of  which  give  the  orig 
Inal,  approximate  Katto-Soda  solution.  In  addition,  solutions 
obtained  numerically  by  digital  computations  are  presented  in 
graphical  and  tabular  form  for  eccentricity  ratios  from  0^0.9 
and  compressible  bearing  parameter,  from  O'*00. 

Design  charts  based  on  these  calculations  are  provided. 
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Symbol 

A 

A' 

B 

B» 

c 

S 

cf 

0L 

D 

h 

k 

L 

M 


N 

P 

q 

R 

s 

T 

u 

v 


Pr 


Meaning 

function  of  11  P* "  defined  by  eq.  68 
function  of  ,,P>U  defined  by  eq.  71 
function  of  "  &  "  defined  by  eq.  65? 
function  of  "  ^ 11  defined  by  eq.  72 
bearing  clearance 
specific  heat  of  fluid 
friction  coefficient  defined  by  eq  .  80 
load  coefficient  defined  by  eq.  78 
diameter  of  shaft 

film  thickness j  also,  fluid  enthalpy 
thermal  conductivity  of  fluid 
length  of  very  long  bearing 

Mach  number  corresponding  to  surface  velocity  of 
shaft;  also,  moment  required  to  turn  shaft 

Prandtl  number  of  fluid  =  cp/«/k 

pressure 

space-dependent  heat  addition  per  umit  volume  per 
unit  time 

radius  of  shaft 

surface  velocity  of  shaft 

temperature 

velocity  of  fluid  in  x-direction;  also,  u  =  7* 

velocity  of  fluid  in  y-direction;  also,  v,  stands 
for  either  sin(k^)  or  cos  (It /S  ) . 

load  on  bearing 


iii 


e  'iV 
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Symbol 


Meaning 

coordinate  in  Cartesian  system;  also,  R© 
coordinate  in  Cartesian  system 
coordinate  in  Cartesian  system 
angle  defined  by  eq.  5*t 

angle  corresponding  to  maximum  film  pressure 

eccentricity,  defined  by  eq.  81 

ff  ag  at  the  point  of  maximum  film  pressure. 

•y  S**-  See  eq.  12. 


angle  around  bearing  periphery,  measured  counter¬ 
clockwise  from  point  of  minimum  film  thickness. 

(ph)niax  See  eq.  i+7. 

fluid  viscosity 

a  Katto-Soda  bearing  parameter,  defined  by  eq.  57. 
also,  kinematic  viscosity  of  fluid. 

mass  density  of  fluid 

l^£cos(y0)  See  eq.  55. 

J  at  the  point  of  maximum  pressure. 

See  eq.  17. 


See  eq.  37. 


dimensionless  pressure  defined  by  eq.  53. 
angular  velocity  of  shaft 
t  Meaning 


(super. ) 
(sub. ) 


denotes  standard  reference  quantity 
zeroth  order,  etc. 
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REFINEMENTS  OF  THE  THEORY  OF  THE  INFINITELY- LONG, 
SELF-ACTING  GAS -LUBRICATED  JOURNAL  BEARING 

Harold  G.  Elrod,  Jr. 

Albert  Burgdorfer 

Introduction: 

Gas-bearing  technology  has  now  a  forty-five  year  his¬ 
tory,  dating  back  to  the  pioneer  work  of  Harrison^ ).  More¬ 
over,  in  this  field  there  has  been  a  recent  spurt  of  activity 
caused  by  the  possibility  of  important  applications.  Never¬ 
theless,  In  the  course  of  progress  towards  useful  devices, 
gaps  In  our  knowledge  have  appeared,  even  in  connection  with 
the  more  elementary  forms  of  lubricating  equipment.  The 
present  paper  endeavors  to  close  a  few  such  gaps  which  have 
appeared  In  connection  with  the  infinitely -long,  self-acting, 
gas- lubricated  journal  bearing.  The  theory  of  this  type  of 
bearing  is  reviewed,  mathematical  techniques  for  numerical 
development  of  the  theory  are  discussed,  and  finally,  im¬ 
proved  design  charts  for  engineering  use  are  presented. 

No  attempt  will  be  made  here  to  review  and  assess  the 
extensive  contributions  of  prior  workers  on  the  subject  of 
this  paper.  A  recent  comprehensive  bibliography of  gas- 
lubrication  research  is  available,  to  which  the  interested 
reader  can  refer.  Two  papers  connected  especially  closely 
with  the  present  research  are  those  of  Ausman(3)  and  of  Katto 
and  Soda(^).  In  each  of  these  papers,  approximate  solutions 
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of  the  gas-bearing  equations  are  given,  together  with  design 
charts  based  on  these  solutions.  Reference  to  these  works 
will  be  made  subsequently  when  comparisons  are  made  with  the 
present  results. 

Basic  Equations  of  Laminar  Lubrication : 

The  general  differential  equations  expressing  the  con¬ 
servation  of  momentum,  mass  and  energy  for  an  arbitrary  New¬ 
tonian  fluid (5)  constitute  the  foundation  for  the  present 
analysis.  Continuity  of  velocity  and  temperature  is  assumed 
at  all  fluid- solid  interfaces.  In  a  paper  presented  in  this 
Symposium,  the  validity  of  these  differential  equations  for  a 
Newtonian  fluid  is  examined  by  Reiner,  and  second-order  cor¬ 
rections  are  proposed.  The  assumed  boundary  conditions  are 
justifiable  in  most  lubrication  applications,  but  Burgdorferi) 
has  shown  that  some  gas  bearings  operate  in  the  slip-flow 
regime.  If  the  molecular  mean  free  path  does  not  exceed 
about  five  per  cent  of  the  radial  clearance  in  the  bearing, 
discontinuity  effects  at  fluid- solid  interfaces  are  not  quan¬ 
titatively  very  significant. 

In  Appendix  I  the  basic  equations  of  laminar  lubrication 
are  derived  from  the  general  conservation  equations  by  means 
of  a  small-parameter  technique  in  which  the  fluid  film  thick- 
ness-to-iength  ratio  (h/L)  serves  as  the  small  parameter.  The 
resulting  equations  are  as  follows  (See  Nomenclature  for  mean¬ 
ing  of  symbols): 
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(1) 

(2) 

(3) 


The  fluid  properties  in  these  aquations  may  be  variable.  Thus 
both  liquids  and  gases  are  included. 

Temperature  Variation  within  Lubricating:  Film: 

The  strict  solution  of  any  lubrication  problem  would  in¬ 
volve  the  coupling  of  the  foregoing  differential  equations  for 
the  fluid  with  the  differential  equations  of  elasticity  and  of 
heat  conduction  appropriate  for  the  shaft  and  bearing.  To  sim¬ 
plify  the  system  of  equations  which  must  be  considered,  we 
shall  here  neglect  both  strain  and  temperature  gradient  within 
the  shaft  and  bearing,  and  assume  that  shaft  and  bearing  are  at 
the  same  temperature.  In  cases  where  this  latter  t-s sumption  is 
Justifiable,  the  Internal  temperature  rise  of  the  film  is  usual' 
ly  negligible,  as  will  now  be  proved. 

To  investigate  the  magnitude  of  the  transverse  temperature 
rise,  we  consider  the  point  of  maximum  film  pressure.  At  this 
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point : 

£$  ~  M  =  ° 

(6) 

In  the 

absence  of  pressure  gradient,  eqs.  1, 

3  and  4  simplify 

to: 

11 

o 

(7) 

(8) 

V*}  <9) 

Let  us  suppose  that  the  surface  at  y  =  0  has  velocity 
components  U,  0,  W  in  the  x,  y  and  z  directions,  respectively, 
and  that  the  surface  at  y  =  h  is  stationary.  The  appropriate 
solutions  of  eqs.  7  and  8  are  then  readily  found  to  be: 


«•  =  V(l-jr) 

ur-  W(>-  T  ) 


Substitution  of  eqs.  10  and  11  into  eq.  9  gives: 

jL  fe  dj  1/  Jjc  i  /Y  dus- 


do) 

(ID 

(12) 

(13) 


4 


THE  FRANKLIN  INSTITUTE  *  Laboratories  for  Research  and  Development 

I-A20L9-10 


Integration  jrields: 

—  -A  U  -h  CO~  +  CoriiSfortlt 

Using  eqs.  1C  and  11  again,  we  gets 


(14) 


Jt  dr 

A*' 


^  +  Conj-fctrit 


where:  32  =  U2  +  W2 

Now  at  the  point  of  maximum  pressure,  oar.  he  re- 
garded  as  a  funotion  of  temperature  only.  Hence  we  define 
the  new  fluid  property: 


(15) 

(16) 


(17) 


£ 


Hi  » 

Because  both  "k"  and  "/l"  are  inherently  positive,  <P  increas¬ 
es  monotonically  with  “T”.  With  the  temperatures  at  y  =  0  and 

V  /  * 

y  =  h  equal,  the  net  change  in  y  across  the  film  is  zero. 

Hence  the  Integral  of  the  right-hand  side  of  sq.  15  must  vanish. 
The  constant  is  thus  determined,  and  the  final  expression  for 

the  distribution  of  5^  is: 


S?±_ 

A 


O-V 


(18) 


^  J  // 

The  maximum  value  of  <p  ,  which  corresponds  to  the  maximum 
value  of  temperature,  occurs  where  ~ 


*L." 


5' 

S' 


(19) 
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Equation  19  is  perfectly  general  for  the  maximum  value 

u 

of  y*  at  the  point  of  maximum  pressure  in  the  film.  Via  have 

t/* 

not  proved  that  this  aquation  gives  the  maximum  value  of  r 
occurring  anywhere  in  the  bearing.  However,  it  does  give  the 
value  of  5^  at  a  very  important  point,  and  this  value  must  be 
at  least  representative  of  values  occurring  elsewhere  in  the 
film. 

An  examination  of  fluid  properties  discloses  (with  no  ob¬ 
served  exceptions)  that  increases  with  temperature  for  both 
liquids  and  gases.  Hence  is  less  in  value  at  the  film 

edges  than  internally.  As  a  result,  we  can  write  the  following 
useful  inequality  for  checking  transverse  temperature  variation 
in  a  lubricating  film. 


Tw  *  *(%), 


(20) 


In  the  case  of  a  perfect  gas  with  constant  specific  heats, 
the  foregoing  inequality  can  also  be  written  as: 


AT 


71 


V  AL  A*  *  M 


(21) 


where  "M"  is  the  surface  Mach  number;  i.  e.,  the  ratio  of  the 
surface  velocity  to  the  velocity  of  sound  in  the  gas. 
Isothermal  Film  Squat ions: 

In  many  applications,  eq.  20  (or  21)  can  be  used  to  es¬ 
tablish  that  the  transverse  temperature  variation  within  the 
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lubricating  film  is  negligible.  In  this  case,  the  temperature 
everywhere  within  the  film  will  be  constant  when  the  bearing 
and  shaft  are  uniformly  at  the  same  temperature.  Hence  all 
fluid  properties  become  functions  of  pressure  only.  Subject 
to  this  condition,  the  differential  equation  for  the  pressure 
in  a  lubricating  film  has  been  derived  many  places.  We  repeat 
the  derivation  here  only  for  the  sake  of  completeness. 

In  an  isothermal  film  all  fluid  properties  are  independent 


of  "y" ,  since  "p11  is  Independent  of  "y". 


become : 


?u.  =• 

dy2-  /*■  Tx 
_  Ji  dfc 

df  2^ 


Thus  eqs.  1  and  3 

(22) 

(23) 


The  boundary  conditions 


u(0)  =  U;  u(h)  =  0 
w(0)  =  Wj  w(h)  =  0 

are  satisfied  by  the  following  solutions  of  eqs.  22  and  23* 


(24) 


*=  +  (25) 

w~=  y/ (/-.£J-h  )  (26> 

The  mass  continuity  equation  (eq.  5)  can  be  integrated 


with  respect  to  "y"  to  give: 

h 

O 


(2?) 
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Equation  27  can  alternatively  be  written  as; 

h  4 

£pfud?  +■ 

o  o 


=  o 


(28) 


because  (a)  the  integration  commutes  with  both  derivatives 

(b)  the  velocities  vanish  at  the  upper  limit,  y  =  h 

(c)  the  fluid  density  is  independent  of  “y" . 

We  can  now  integrate  the  velocities  given  by  eqs.  25  and  26, 
as  required  in  aq.  28,  to  getr 


Or; 


>(v*f£*w%P  "" 


There  is  a  useful  exact  analogy  which  should  here  be 
pointed  out  between  the  pressure  in  a  lubricating  film  and 
the  temperature  in  a  constant  density,  moving  fluid.  To  il¬ 
lustrate  this  fact,  let  us  rearrange  eq.  30  in  the  following 
form; 


In  this  equation,  "p"  plays  the  role  of  temperature,  while 
"U”  and  "W"  retain  their  meaning  as  velocities.  The  quantity 
fa  ^4)^Is  intrinsically  positive,  because  the  film  thickness 
is  certainly  positive,  and  *  ^he  "Newtonian  velocity  of 
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sound,0  Is  positive  for  any  mechanically  stable  substance. 
Hence  we  let  6$ play  the  role  of  volumetric  specific  heat. 

The  quantity  is  inherently  positive,  so  we  let  is  serve 

in  the  analogy  as  thermal  conductivity.  Finally,  the  last 
group  of  terms,  being  the  product  of  "p  u  with  a  space-depen- 


dent  function  £ p  ^ + JV )  *  we  take  t0  be  a  temperature  and 
space-dependent  Internal  source  of  heat.  Thus  the  equation  in 
convective  heat  transfer  analogous  to  eq.  31  is* 

+ f(r>^ 


(32) 


The  terms  in  this  last  equation  have  been  defined  by  the  pre¬ 
ceding  discussion.  This  analogy  enables  one  to  exploit  intui¬ 
tion  and  solutions  acquired  in  another,  well-established  field. 

Differential  Equation  and  Mass  Content  of  Infinitely -Long 
Self-Acting  Journal  Bearing: 

Figure  1  shows  schematically  a  long,  self-acting  journal 

bearing.  Here  "x"  denotes  peripheral  distance  measured  in  a 

counter-clockwise  direction,  and  "z"  denotes  axial  distance. 

The  use  of  "x”  to  denote  a  curvilinear,  rather  than  recti- 

(7) 

linear  distance,  has  been  examined  by  Elrod  for  the  case 
of  fluids  with  constant  properties,  and  the  effect  has  been 
found  to  be  exceedingly  slight.  There  Is  no  reason  to  suppose 
that  otherwise  will  be  the  case  for  fluids  with  variable  proper¬ 
ties. 


When  a  journal  bearing  is  perfectly  aligned,  with  both  ends 
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exposed  to  the  same  ambient  pressure,  there  is  a  simple  conserva¬ 
tion  condition  which  can  be  derived  from  eq.  30.  For  this  situ¬ 
ation:  h  =  h(x)  and  W  =  0  (33) 

and  eq.  30  reduces  to: 


<■3 
7. 

Now  all  flow  quantities  are  cyclic  in  nxM;  i.  e.,  their  net 
change  circumferentially  around  the  shaft  is  zero.  Hence,  if 
we  integrate  eq.  34  around  the  shaft,  the  right-hand  side  re¬ 
duces  to  zero,  and  we  are  left  with: 


(34) 


Or: 


H  dx  dx  - 0  (35) 

/if  J&.  dx  —  A Cl  Constant  (36) 

Since  both  "f"  and  "/ t"  are  functions  of  "p"  only,  we 


can  define  a  new  function:  pX^  — '  f 
in  terms  of  which  eq.  36  becomes: 


(  37) 


(38) 


*  *•  ^ b  ~  f  %  t-  &  (39) 

The  cyclic  integral  in  eq.  39  must  have  the  same  value  at 
z  =  0  and  z  =  L.  Hence  the  constant  11  A*'  must  be  zero.  The 
following  conservation  equation  then  results: 
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(40) 


where  corresponds  to  ambient  pressure,  and  may  be  taken 

as  zero,  if  desired.  Sines,  very  closely, 

h  =  c  (1  +  e  cos  0  )  (41) 

we  haves 

^/?3‘dp  —  cK  *§(/-{- 6  CM  e)  d&  ~  ZirKc  (i+ze  )  (42) 

$h3p6<tx  =  ZirRcCt+ieV^C^  (43) 


This  last  equation  applies  to  any  Newtonian  fluid,  in  the  ab¬ 
sence  of  cavitation. 

If  the  lubrioant  is  a  perfect  gas,  then  its  viscosity  is 
independent  of  pressure,  and,  at  constant  temperature,  its  den¬ 
sity  is  proportional  to  the  pressure.  Hence  we  can  take  /\- 
equal  to  p2.  That  isj 

j>  tf/fdtx  =■  coolant  =  C  (t  +  )  (£  (44) 

Since  this  condition  applies  to  a  gas  journal  bearing,  how¬ 
ever  long,  it  applies  to  the  "infinitely  long”  bearing. 

That  is,  it  applies  to  the  central  region  of  a  long  bearing, 
far  from  either  end,  where  the  pressure  is  essentially  a  func¬ 
tion  of  "x"  only.  It  is  only  in  terms  of  this  region  that  we 
may  think  of  the  "infinitely  long  bearing",  since  otherwise  we 
can  establish  no  connection  with  an  ambient  pressure,  and  hence 
no  connection  with  any  physically  real  problem. 
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Equation  43,  with  its  various  specializations,  such  as 
eq.  44,  furnishes  a  condition  which  sets  the  general  pressure 
level,  and  hence  the  mass  content,  of  any  continuous-film,  in¬ 
finite  bearing.  It  is  an  important  relation  in  the  theory  of 
the  infinite  gas  bearing,  and  has  not,  so  far  as  the  authors 
are  aware,  been  previously  employed.  In  the  case  of  non-cavi- 
tating  films  of  incompressible  fluid,  the  mass  content  of  the 
bearing  is,  of  course,  known.  Also,  if  the  viscosity  is  con¬ 
stant,  the  load-carrying  capacity  of  the  bearing  can  be  com¬ 
puted  without  knowledge  of  the  general  pressure  level.  How¬ 
ever,  eq,  43  still  applies,  and  might  be  useful  for  predicting 
the  onset  of  cavitation,  or  for  refined  treatments  with  variable 
viscosity. 

The  differential  equation  for  the  infinitely-long  bearing 
is  obtained  from  eq.  34  by  eliminating  z-derivatives.  It  is: 

d  <&  =  6  U  (45) 

dx  A  dx  dx 

Integrating  once,  we  obtain  the  following  first-order  non¬ 
linear  differential  equation  whose  solution  involves  specifi¬ 
cation  of  two  constants,  one  of  which  Is  shown  explicitly. 

ph  dfi  —  1/  flh  +  Cons  fay  (46) 

yfc  dx  ' 

The  two  constants  are  determined  by  requiring  (a)  that  the 
solution  be  cyclic  in  "x"  and  (b)  that  it  satisfy  the  condi¬ 
tion  given  by  eq.  43. 
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Series  Solution  for  the  Pressure  in  Infinitely- Long ,  Self-Act¬ 
ing  'Gas-Lubricated  Journal  Bearing 

For  gas-lubricated  bearings,  eq.  *f6  can  be  readily  put  in 


where  6  =  x/R  and  "  k  11  is  the  value  of  the  product  nph"  at  a 
stationary  point  for  the  pressure.  When  the  pressure  in  the 
film  varies  and  the  film  thickness  is  given  by  eq.  ^1,  there 
are  two  such  stationary  points,  one  corresponding  to  the  maximum 
pressure  and  one  corresponding  to  the  minimum  pressure  "P2". 


Since  "a  11  Is  limited  to  a  single  value,  we  must  have; 


•45 

< 

II 

■4T* 

< 

(48) 

Or: 

(49) 

Now  the  ratio  of  h^/h2  is  certainly  less  than  the  maximum 
value  of  film  thickness  c(f  +  e)  divided  by  the  minimum  value 
of  film  thickness  c(f  -e).  Hence: 


me  xorm: 


—  =s  L±  &  <5o) 

I  £ 

Prom  eq.  bk-  it  is  easily  seen  that  p^  is  greater  than  pa  and 
that  p0  is  less  than  p  .  Then  from  eq.  50  we  get: 

Cm  Cl 


(5D 

(52) 

It  is  evident  from  eq.  52  that  the  load-carrying  capacity  of 
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a  gas-lubricated  journal  bearing  has  a  limit  which  cannot  be 
exceeded  through  changes  in  gas  or  in  speed,  but  only  through 
changes  in  allowable  eccentricity  or  in  ambient  pressure. 

Closed-form  solutions  of  eq,  47  do  not  exist,  but  approxi¬ 
mate  analytical  solutions  can  be  very  helpful.  Both  Ausman^) 

(41 

and  Katto  and  Soda  have  developed  such  solutions.  In  the 
analysis  which  follows,  a  series  expansion  for  the  film  pres¬ 
sure  will  be  obtained,  the  first  terms  of  which  give  the  Katto 
and  Soda  solution,  with  succeeding  terms  giving  refinements 
thereof.  Let  us  begin  by  introducing  the  Katto  and  Soda  vari¬ 
ables.  They  arej 


qp  5?  fc£:(i~e) 

i  /c 

(53) 

,  e  +  C44,  0 

6*6/3  -  - - ^ 

W 

In  terms  of  these  new  variables,  the 

differential  equation 

corresponding 

to  eq.  4 7  is: 

/  rr,  _  s. ' 

d/»  ~  7  TU  rj 

7  (55) 

where : 

—  /  —  B  C04> f 

(56) 

and: 

(57) 

6ft  l/R 

Examining  eq.  55,  we  note  that  when  e  =  0,  the  only  cyclic 
solution  is  'f  =  1  =  t*  .  On  the  other  hand,  when  the  parameter 
O  ,  we  know  the  pressure  variation  must  nevertheless  be 
constrained  by  relations  50-52.  Then  if  the  pressure  derivatives 
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are  also  to  be  bounded,  we  observe  from  eq.  55  that  must 

again  approach  "TM.  In  view  of  these  facts,  we  hypothesise 
u  =  T  as  a  small  quantity  and  attempt  an  expansion  In  the 


formi  u  =  e  ut(/z)  4-  €?uz(j3)  4-  g3u3(/3)  +  ' 

In  terms  of  "u”  the  differential  equation  becomes: 

+  Zr  +  2  U.  p-  +24*  p'  —  -f  *ru 
dfi  */$  dj3  * 

Now  in  terms  of  the  series  expansion  in  eq.  58, 


(58) 


(59) 


u. 


<x>  ext 

l+J  » 


-  -Z* 

v 


(60) 


J  —  ° 


Here  the  convention  is  adopted  that  =  0  for  n  £0. 
Using  the  expansions  from  eqs.  58  and  60  in  eq.  59,  we  get: 


+J?  (61) 


n=0 


+  *^*±2^,  +z£<?XnW= 


»  SO 


nzs-o 


n=0 


n-o 


where, o  ;  %  ^  =/W/3 ;  ^ (fi)= 

%„  (fl)  =  o,  «23  (62) 

When  the  coefficients  of  the  various  powers  of  " €  "  are 
equated  to  zero  in  eq.  61,  the  following  sequence  of  linear 
differential  equations  results. 

A  K  -  &&  u»-i)  /eetiA  ~  ^ 


Oo 


-X  W -it  UJ  “»-j 

J-& 

16 


(63) 


THE  FRANKLIN  INSTITUTE  *  Laboratories  for  Research  and  Development 

I-A20L9-10 


The  cyclic  solution  for  (un  -  ^n-i003/9  )  is  given  by: 

<**/««-«  +  «*> 


where  the  constant  of  integration  is  to  be  taken  as  zero  to 
avoid  an  exponentially-increasing,  non-periodic  solution. 

The  products  i^u^j  can  all  be  reduced  trigonometrically 
to  a  sum  of  terms  of  the  form  cos(k^)  or  sin(k^),  where  ’’k" 
is  an  integer.  With  vk(/3)  standing  for  either  of  these  func¬ 
tions,  the  following  indefinite  integral  formula  is  easily 
derived: 

"  e  *> ' <4  &  df>  =  aT^/s  P  (6?) 

With  the  aid  of  this  formula  we  find  from  eq.  64: 


u,  — 


/V-^: 


^  ~h  ^  ~j 


U  — 


-  £ 


/  * 


$‘}d/* 


(65) 

(66) 


By  squaring  eq.  65  we  find  the  terms  in  the  curly  brackets 
of  the  integral  in  eq.  66  to  be: 


+  i 


A  A****  4-  /0*i< 


(67) 


where : 


A  - 
6  = 


J_  ±  A*  0-^^ 

Z  (f +**’)*' 
l>3 

(/+*>*-)*- 


(68) 

(69) 


17 


THE  FRANKLIN  INSTITUTE  <■  Laboratories  for  Research  and  Development 

I-A20L9-10 


Applying  eq.  66  again,  we  get: 


—  i  +  [(A  ~  Z.Z  Q)  -f-  {5-4  Z  AA)yC&& 


a  -  zp’b = a'  =  -±  r/y;^r; 

B  +  Z^A  =  B'~  -t>  [  /+-2Z4-  ~j 

L  o+pv*-J 


(70) 


(7D 

(72) 


Substituting  the  expressions  for  u^  and  U2  in  series  58,  we 
obtain  the  following  expression  for  valid  to  0(e2) . 

ft  =  /  —  £■  /!&4sy6  ~f~  £  (/  -f-  €  •+■  tZjiCo&A )  ~t~ 

(a'a^~Z/>  -h  b's^Zjs)  +■ .  (73) 

To  0(e)  eq.  73  gives  precisely  the  Katto-Soda  solution; 


i.  e. 


ft  ~  / 


(74) 


Typical  circumferential  profiles  of  the  dimensionless 
pressure  "  are  shown  in  Pigs.  2  and  3*  as  computed  from 
eq.  73*  These  profiles  are  compared  with  others  obtained  by 
an  accurate  numerical  method,  presently  to  be  described.  It 
will  be  observed  that  the  series  solution  is  most  accurate 


for  small  "  i>  11  and  small  "e".  This  characteristic  was  anti¬ 
cipated  in  the  formulation  of  series  58.  It  will  also  be  ob¬ 
served  that  the  series  solution  is  scarcely  of  adequate  accur¬ 
acy  for  the  operating  condition  defined  by  Z  =  1,  e  -  0.8  . 
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Digital  Computer  Solutions : 

In  view  of  the  difficulty  in  finding  approximate  solutions 
of  the  one-dimensional  gas  bearing  equations  with  sufficient 
accuracy,  it  was  decided  to  attack  the  problem  directly  by 
using  standard  numerical  techniques  to  solve  eq,  55.  The 
method  employed  was  the  interative  scheme  of  Cllppinger  and 

o 

Dims dale  .  The  angle  Vo"  corresponding  to  a  maximum  value 
of  ^  was  estimated  from  the  ICatto  and  Soda  solution. 

Then  the  solution  was  started  at  this  angle  and  continued 
in  a  clockwise  direction  with  an  angular  increment  of  7r/30 . 
Three  "rotations"  of  the  computer  usually  sufficed  to  produce 
a  solution  periodic  to  within  six  significant  figures. 
Dimensionless  pressure  profiles  were  obtained  for  five  differ¬ 
ent  eccentricities  (e  =  0.2,  0.4,  0.6,  0.8,  0.9)  with  each 
of  four  different  values  of  the  speed  parameter  "jSn  (  ~  1,00, 

0.5,  0.25,  0.125).  Values  of  "  are  listed  in  Tables  2-5. 
They  are  believed  to  be  correct  to  within  one  digit  in  the 
last  figure  tabulated. 

The  accuracy  of  both  the  series  and  digital  computer  com¬ 
putations  is  confirmed  by  the  close  agreement  between  the  re¬ 
sults  of  both  for  small  "  and  "£  ",  as  shown  in  Pigs.  2 
and  3,  Actually,  the  derivation  of  the  series  was  purely 
formal,  no  proof  of  convergence  being  offered.  As  a  check  on 
convergence,  the  following  case  was  computed  by  both  methods: 
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l?  =  1.0  ;  6  =  0.05 

The  results  are  compared  belov;  in  the  neighborhood  of  the 
pressure  extremes*  It  will  be  observed  that  the  inclusion 

of  terms  beyond  the  original  Katto-Soda  solution  does  yield 
a  significant  improvement  In  precision,  tending  to  confirm 

the  assumption  of  convergence  of  the  series. 


TABLE  1  : 


COMPARI SOU  OP  DIGITAL  COMPUTER  AND  SERIES  SOLUTIONS 


FOR  THE  CASS  ^  =  1.0; 

6  -  0.0 5 

Values  of  " 

/8  (rad) 

Eq.  74 

^Z3L 

Computer 

2.40855 

2.35619 

2.30383 

1.03531 

1.03536 

1.03531 

1.03583 

1.03586 

1.03578 

1.03580 

1.03583 

1.03575 

5.55015 

5.49779 

5.44543 

0.96469 

0.96465 

0.96469 

0 196522 
0.96514 
0.96517 

0.96525 

0.96518 

0.96520 

20 


NUMERICAL 


COMPARISON^  OF  IMPROVED,  K-S  SOLUTION 
WITH  EXACT  NUMERICAL  SOLUTION  V‘1.00 


COMPARISON  OF  IMPROVED  K-S  SOLUTION 
WITH  EXACT  NUMERICAL  SOLUTION  V*  0.125 
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PERIODIC  SOLUTIONS  OF  THE  KATTO-SODA  DIFFERENTIAL  EQUATION 

rr  r  7- 

=  1.00 


e 

0.20 

0-40 

0.60 

0.80 

0.90 

cosx-9" 

-0.707 

-0.629 

-0.629 

-0.629 

-0.629 

iff  AS) 

i* 

-t 

*4* 

0 

1.148 

1.300 

1.459 

1.621 

1.703 

l 

1.144 

1.285 

1.435 

1.589 

1.667 

2 

1.134 

1.261 

1.399 

1.541 

1.613 

3 

1.119 

1.228 

1.352 

1.481 

1.547 

4 

1.100 

1.189 

1.298 

1.413 

1.473 

5 

1.078 

1.146 

1.239 

1.343 

1.399 

6 

1.053 

1.099 

1.178 

1.274 

1.329 

7 

1.026 

1.052 

1.118 

1.212 

1.269 

8 

0.999 

1.004 

1.061 

1.158 

1.221 

9 

0.972 

0.958 

1.008 

1.113 

1.186 

10 

0.946 

0.914 

0.957 

1.075 

1.162 

11 

0.928 

0.873 

0.910 

1.041 

1.142 

12 

0.902 

0.836 

0.863 

1.006 

1.121 

13 

0.886 

0.805 

0.820 

0.966 

1.092 

it 

0.874 

0.783 

0.780 

0.922 

1.053 

15 

0.869 

0.771 

0.751 

0.875 

1.004 

16 

0.870 

0.774 

0.738 

O.836 

0.953 

17 

0.878 

0.791 

0.748 

0.816 

0.915 

18 

0.892 

0.823 

0.784 

0.829 

0.906 

19 

0.912 

0.868 

0.842 

O.878 

0.937 

20 

0.937 

0.922 

0.918 

0.959 

1.008 

21 

0.966 

0.982 

1.004 

1.059 

1.106 

22 

0.996 

1.043 

1.094 

1.169 

1.220 

23 

1.027 

1.104 

1.183 

1.279 

1*337 

24 

1.056 

1.160 

1.126 

1.381 

1.447 

25 

1.083 

1.209 

1.335 

1.471 

1.544 

26 

1.107 

1.249 

1.392 

1.543 

1.622 

27 

1.125 

1.279 

1.434 

-1.595 

1.679 

28 

1.139 

1.298 

1.459 

1.625 

1.711 

29 

1.146 

1.305 

1.467 

1.634 

1.719 

Note;  ^>o  is  in  the  second  quadrant. 
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PERIODIC  SOLUTIONS  OP  THE  KATTO-SODA  DIFFERENTIAL  EQUATION 

£3*-  -r/7  rl  . 
d/$  ~  »(_  'V'J  * 

^  =  0.50 


/  -  £ 


<5  — ► 

cos  fo-*- 

) 

0*20 

-0.843 

nj. 

0.40 

-0.843 

0 

1.179 

1.360 

1 

1.171 

1.343 

1.314 

1.272 

1.220 

2 

l 

1.156 

1.135 

1.108 

5 

1.076 

1.160 

6 

1.042 

1.094 

7 

1.006 

1.026 

8 

0.970 

0.958 

9 

0.935 

0.891 

10 

0.902 

0.828 

11 

0.874 

0.772 

12 

0.851 

0.724 

13 

0.835 

O.687 

14 

0.825 

0.664 

15 

0.824 

0.657 

16 

0.831 

0.669 

3.7 

0.847 

0.698 

18 

0.869 

0.744 

19 

0.898 

0.803 

20 

0.931 

0.871 

21 

0.968 

0.945 

22 

1,006 

1.021 

23 

1.043 

1.096 

24 

1.079 

1.167 

25 

1.111 

1.229 

26 

1.138 

1.282 

27 

1.159 

1.323 

28 

1.173 

1.350 

29 

1.180 

1.362 

0.60 

0.80 

0.90 

-0.843 

-0.843 

-0.843 

V 

V 

'll* 

1.542 

1.724 

1.815 

1.317 

1.472 

1.690 

1.777 

1.631 

1.7H 

1.410 

1.551 

1.621 

1.335 

1.453 

1.512 

1.239 

1.343 

1.391 

1.156 

1.227 

1.265 

1.061 

1.111 

1.142 

0.967 

1.002 

I.029 

0.877 

0.903 

0.933 

0.793 

O.818 

0.858 

0,716 

0.745 

0.803 

0.648 

O.683 

0.762 

0.589 

0.624 

0.725 

0.544 

0.566 

0.681 

0.521 

0.511 

0.624 

0.526 

0.475 

0.563 

0.563 

0.484 

0.528 

0.629 

0.547 

0.553 

0.718 

0.652 

0.643 

0.821 

0,784 

0.777 

O.932 

1.046 

0.928 

0.933 

1.077 

1.096 

1.157 

1.223 

1.258 

1.460 

1.353 

1.359 

1.479 

1.409 

1.544 

1.430 

1.580 

1.656 

1.489 

1.657 

1.742 

1.528 

1.708 

1.798 

1.546 

1.730 

1.823 

Note:  A  is  in  the  second  quadrant. 
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TABLE  4 


PERIODIC  SOLUTIONS  OF  THE  KATTO-SODA  DIFFERENTIAL  EQUATION 
J  'll) 

*  -  -1-  '  ' 


_  rr.  _rl  . 

l>=  0.25 


£  — *■ 
cos  /Q0  —* 

£(Ar/S 

0 

1 

2 

l 

5 

6 


9 

10 

u 

12 


15 

16 

17 

18 

19 

20 
21 
22 

3 

25 

26 
2? 
28 
29 


0.20 

-0.940 

V' 

1.193 

1.185 

1.168 

1.145 

1.115 
1.080 
1.042 
1.002 
0.962 
0.924 
0.888 
0.858 
0.833 
0.816 
0.808 
0.807 
0.816 
0.832 
0.856 

O.887 

0.923 

0.962 

1.002 

1.042 

1.081 

1.116 
1.145 
1.169 
1.185 
1.193 


0.4o 

-0.940 

1.386 

1.369 

1.337 

1.290 

1.232 

1.163 

1.088 

1.009 

0.930 

0.854 

0.783 

0.721 

0.672 

0.636 

0.617 

0.615 

0.632 

0.666 

0.715 

0.777 

0.849 

0.927 

1.008 

1.088 

1.165 

1.234 

1.293 

1.339 

1.371 

1.387 


0.60 

-0.940 

1.580 

1.554 

1.506 

1.437 

1.349 

1.248 

1.137 

1.021 

0.905 

0.793 

0.690 

0.598 

0.522 

0.465 

0.432 

0.426 

0.450 

0.501 

0.576 

0.670 

0.779 

0.896 

1.017 

1.137 

1.251 

i:$ 

1.510 

1.557 

1.581 


0.80 

-0.940 

1.773 

1.740 

1.675 

i:83 

1.335 

1.190 

km 

0.747 

0.618 

0.506 

0.412 

0.335 

0.279 

0.253 

0.275 

0.342 

O.443 

0.569 

0.713 

0.868 

1.029 

1.188 

1.338 

1.475 

1.591 

1.682 

1.744 

1.775 


0.90 

-0.956 

1.874 

1.845 

1.781 

1.684 

1.562 

1.417 

1.258 

1.092 

0.926 

0.768 

0.628 

0.509 

0.417 

0.347 

0.289 

0.235 

0.209 

0.250 

0.350 

0.484 

0.640 

0.811 

0.990 

1.170 

1.342 

1.499 

1.635 

1.745 

1.823 

1.867 


Note:  /3a  is  in  the  second  quadrant. 
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PERIODIC  SOLUTIONS  OF  THE  KATTO-SODA  DIFFERENTIAL  EQUATION 

-%['-£] ;  Tm  !-***«/> 

l>  =  0.125 


€  — *- 

cos a  — 

IF  (Art) 

0.20 

-0.992 

0.4-0 

-0.992 

0.60 

-0.992 

_  ^  ... 

O.80 

-0.992 

.  V 

0.90 

-0.992 

_  tP  ... 

0 

1.198 

1.397 

1.595 

1.794 

1 

1.194 

1.388 

1.582 

1.777 

2 

1.181 

1.363 

1.544 

1.726 

1.817 

3 

1.161 

1.322 

1.483 

1.644 

1.724 

4 

5 

1.133 

1.266 

1.400 

1.534 

1.601 

1.100 

1.200 

1.300 

1,402 

1.452 

6 

1.062 

1.125 

1.188 

1.252 

1.284 

7 

1.021 

1.044 

1.067 

1.092 

1.105 

8 

0-980 

0.961 

0.944 

0.929 

0.922 

9 

0.939 

0.880 

0.823 

0.769 

0.743 

10 

0.901 

0.804 

0.710 

0,619 

0.576 

11 

0.868 

0.737 

0.608 

0.486 

0.4-28 

12 

O.84o 

0.681 

0.-524 

0.374 

0*306 

13 

0.819 

0.639 

0.460 

0.287 

0.211 

14 

0.806 

0.612 

0.419 

0.230 

0.146 

15 

0.802 

0.603 

0.405 

0.208 

0.114 

16 

0.806 

0.612 

0.418 

0.225 

0.129 

17 

0.819 

0.638 

0.457 

0.278 

0.190 

18 

O.84o 

0.680 

0.521 

0.364 

0.287 

19 

0.868 

0.736 

0.605 

0.477 

0.414 

20 

0.901 

0.803 

0.707 

0.612 

0.566 

21 

0.939 

0.880 

0.821 

0.764 

0.736 

22 

0.980 

0.961 

0.943 

0.926 

0.918 

23 

1.021 

1.044 

1.067 

1.091 

1.104 

24 

1.062 

1.125 

1.188 

1.252 

1.284 

25 

1.100 

1.200 

1.301 

1.401 

1.402 

1.453 

26 

1.133 

1.267 

1.535 

1.645 

1.602 

27 

1.161 

1.322 

1.483 

1.726 

28 

1.181 

1.363 

1.545 

1.727 

1.818 

29 

1.194- 

1.388 

1.583 

1-777 

1.874 

Note:  A  is  in  the  second  quadrant. 
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In  order  that  the  computed  data  on  pressure  distributions 
my  be  useful,  it  is  necessary  to  relate  the  internal  pressures 
of  the  infinite  bearing  to  some  ambient  pressure  by  means  of 
the  mass  content  rule.  When  a  change  is  made  from  convention- 
al  variables  to  those  employed  by  Katto  and  Soda,  the  mass 
content  rule  (eq.  44)  gives: 


at  & 

=  «**('- **) 

(75) 

=  zrc*(/+ £<*)/}* 

(76) 

Or:  = 

Ac  < 

j  2r0+%  £*) 

±\j/S 

(77) 

For  the  data  tabulated  in  Tables  2-5  the  required  integra¬ 
tion  in  eq.  77  "was  performed  using  Simpson* s  rule  with  an 
angular  interval  of  IT/ 15. 

With  k/ Pqc  known,  the  pressure  ratio  p/pQ  can  readily 

m  a 

be  found  in  terms  of  operating  variables.  Thus: 


p/b  =  *-  jk 

r/F°-  /tp  /--  ■ 


(78) 


/c 


G/tlTR 


C79) 


To  prepare  Design  Charts  in  terms  of  conventional  bear¬ 
ing  parameters,  other  integrations  of  the  -curves  are 
required.  For  these  integrations  Simpson's  rule  was  again 
used  with  an  angular  interval  of  fT/15* 
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Design  Charts i 

The  steady-state  performance  of  an  infinitely-long,  self¬ 
acting,  gas-lubricated  journal  bearing  can  be  completely  ex¬ 
pressed  in  terms  of  certain  well-known  bearing  parameters. 

In  order  that  the  subject  matter  of  this  section  may  be  self- 
contained,  these  parameters  will  here  be  defined  in  terms  of 
physical  quantities  which  are  easily  interpreted.  These  par¬ 
ameters  can  be  readily  found  from  the  Tp  -curves,  with  the 

use  of  eqs.  77-79* 

Load  Parameter,  CL 

This  parameter  represents  the  ratio  of  the  load,  V/,  which 
is  actually  supported  by  the  bearing  to  the  force  which  would 
be  exerted  by  ambient  pressure  acting  over  the  projected  area 
of  the  shaft;  i.  e., 


(78) 


Attitude  Angle, £ 

The  attitude  angle  for  a  shaft  rotating  counter-clock¬ 
wise,  as  shown  in  Fig.  1,  is  the  angle  between  the  load  vec¬ 
tor,  W,  and  the  radius  vector  to  the  point  of  minimum  film 
thickness;  i.  e., 

f  m  Z7T-  (79) 

Friction  Coefficient, 

The  friction  parameter  is  the  ratio  of  the  torque,  M, 
actually  required  to  rotate  the  shaft  in  its  eccentric 
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position  under  load  W,  to  the  torque,  MQ,  which  would  be 
required  to  rotate  the  shaft  if  it  were  concentric  with  the 
bearing;  i.  e., 


9* 


M 


(80) 


Eccentricity  Ratio, 6 

This  parameter  measures  the  lack  of  concentricity  of 
the  shaft  with  the  bearing,  and  is  the  ratio  of  the  difference 
between  maximum  and  minimum  film  thicknesses  (h^y  -  ^^5  to 
the  difference  between  bearing  and  shaft  diameters  (2c);  i.e., 

A  _  L 

£  —  max _ ntr>itt  (81) 

Zc- 

Compressibility  Bearing  Parameter, A. 

This  parameter  is  3/n*  (practically  unity)  times  the 
ratio  of  the  torque  for  concentric  operation,  MQ,  to  the 
torque  that  would  be  produced  by  ambient  pressure  acting 
on  the  average  cross-sectional  film  area,  cL;  i.  e,, 


Table  6  summarizes  the  inter-relationships  of  the  fore¬ 
going  parameters,  as  found  by  digital  computation  for  speci¬ 
fic,  discretely-spaced  operating  conditions.  Design  charts 
based  on  these  data  (suitably  extrapolated  by  the  formulas 
of  Ausman  and  Katto-Soda.)  are  presented  in  Figs.  4,  5  and 
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Table  6 

SUMMARY  OF  DIGITAL  COMPUTER  COMPUTATIONS 


£ 

Cl 

% 

4 

1.0 

0.2 

0.9622 

0.2207 

1.0518 

45.75 

0.4 

0*8209 

0.4156 

1.2346 

47.97 

0.6 

o.54i8 

0.4518 

1.6536 

52.78 

0.8 

0,2052 

0.3039 

2.6722 

54.53 

0.9 

0.0715 

0.1708 

4.0505 

52.88 

0.5 

0.2 

1.979 

0.2914 

1.0331 

26.93 

0.4 

1.859 

0.6367 

1.1528 

28.10 

0.6 

1.460 

0.9967 

1.4501 

30.73 

0.8 

0.6348 

0.9537 

2.3268 

35.05 

0.9 

0.2221 

0.5933 

3.6303 

32.49 

0.25 

0.2 

4.012 

0.3214 

1.0242 

14.03 

0.4 

3.976 

0.7534 

1.1090 

14.38 

0.6 

3.643 

1,440 

I.3090 

15.10 

0.8 

2.379 

2.447 

1.8987 

17.23 

0.9 

1.017 

2.298 

2.9412 

19.80 

0.125 

0.2 

8.058 

0.3302 

1.0214 

7.13 

O.k 

8.110 

0.7869 

1.0956 

7.16 

0.6 

7.781 

1.581 

1.2648 

7.31 

0.8 

0.9 

6.256 

4.246 

3.379 

5.461 

1.7220 

2.4733 

7.66 

8.27 
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6.  With  their  use  the  designer  can  quickly  arrive  at  an 

approximate  bearing  design.  Approximate  corrections  for  the 

differences  in  behavior  between  bearings  of  finite  and  of 

9 

infinite  length  are  given  by  Ausman  in  a  paper  presented 
in  this  Symposium,  ^ven  when  designed  according  to  this 
method,  the  bearing  may  not  perform  satisfactorily  because 
the  predicted  mode  of  steady- state  operation  is  unstable. 

The  theory  of  this  "whirl"  phenomenon  is  not  yet  fully  under¬ 
stood,  so  that  each  bearing  design  should  be  "proven"  by  ex¬ 
periment.  Further  design  considerations  are  given  in  ref.  9. 

Several  features  of  the  curves  in  Fig.  b  are  worthy  of 
comment.  First,  at  high  values  of  the  compressibility  par¬ 
ameter,  A,  the  load  coefficient,  C^,  becomes  dependent  on 
the  eccentricity  ratio, 6  ,  only.  (This  fact  might  have  been 
anticipated  from  the  inequality  in  eq,  52.)  Thus  an  increase 
of  the  rotational  speed  of  a  bearing  may  sometimes  result  in 
a  negligible  increase  of  load-carrying  capacity.  Also,  if 
the  speed,  ambient  pressure,  overall  dimensions  and  minimum 
film  thickness  of  a  bearing  be  maintained  constant,  an  increase 
of  load-carrying  capacity  may  sometimes  be  achieved  through 
an  increase  of  clearance,  cl  Second,  at  small  values  of  the 
eccentricity  ratio  the  attitude  angle  becomes  dependent  on 
the  compressibility  parameter  only.  Thus  the  compressibility 
parameter  sets  the  phase,  so  tc  speak,  of  the  pressure  distrib- 
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ution  within  the  hearing  relative  to  the  film- thickness  dis¬ 
tribution  in  a  mariner  somewhat  similar  to  the  way  the  free- 
stream  Mach  number  affects  the  pressure  distribution  on  a 

10 

wave-shaped  wall.  Third,  at  very  low  values  of  the  com¬ 
pressibility  parameter,  the  behavior  of  the  complete,  liquid- 
filled,  non- cavita ting  journal  bearing  is  approached.  This 
region  is  crowded  into  the  lower  left-hand  corner  of  Fig.  4. 

With  respect  to  Fig.  f>,  which  depicts  the  characteris¬ 
tics  of  the  friction  coefficient,  C^,  it  should  be  observed 
that  for  small  values  of  the  eccentricity  ratio  the  coeffic¬ 
ient  is  nearly  that  for  concentric  operation.  At  A  =  0,  the 
"incompressible1*  coefficients  are  shorn,  as  constructed 
from  the  data  in  ref.  11.  The  important  features  of  the 
curves  of  attitude  angle  in  Fig*  6  have  already  been  discus¬ 
sed  in  connection  with  Fig.  4. 
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APPENDIX  I 


Derivation  of  the  Equations  of  Laminar  Lubrication: 

In  this  appendix  the  basic  equations  of  laminar  lubrica¬ 
tion  are  derived  from  the  over-all  conservation  equations  of 
momentum,  energy  and  mass  for  a  single-phase  Newtonian  fluid. 
The  geometry  treated  is  that  for  the  infinite  slider  bearing. 
Generalisation  of  the  derivation  to  include  finite  slider 
bearings  is  obvious,  and  generalization  to  include  curved 
surfaces  can  no  doubt  be  accomplished  with  general  tensor 
techniques  similar  to  those  employed  in  ref.  7. 

The  physical  situation  to  be  analysed  is  shown  in 
Fig.  7.  The  upper  surface  of  the  bearing  is  stationary. 


Schematic  Diagram  of  Slipper  Bearing 


The  lower  surface  moves  in  the  x-direction  with  velocity  U. 

The  fluid  pressures  at  the  entrance  and  exit  of  the  lubricating 
film  have  the  same  value,  p*.  The  bearing  surfaces  are  taken 
to  be  at  the  same  uniform  temperature,  T*. 

Basic  Equations ; 

The  eauatlons  governing  the  two-dimensional  motion  of 
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a  Newtonian  fluid  are  as  follows  (See,  for  example,  ref.  5) ’ 
Conservation  of  Linear  Momentum: 


Conservation  of  Mass : 

3J£^=  O  (A3) 


Conservation  of  Energy: 


The  boundary  conditions  to  be  satisfied  are  as  follows: 


u  =  U,  v  =  0,  T  =  T*;  y  =  0 

u  =  0,  v  =  0,  T  =  T*;  y  =  h(x)  (A?) 

p  =  p*;  x  =  0  or  L 

Formulation  of  Dimension- free  Equations: 

In  order  to  determine  the  relative  magnitudes  of  the 
terms  in  equations  l-*f  in  the  limiting  case  of  a  very  thin 
fluid  film,  it  is  convenient  to  convert  all  terms  to  appropriate 
dimensionless  variables.  First  we  define  the  new  independent 
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variables:  x  =  x/L  and  f  -  y/h(x)  (A6) 

Then,  as  In  ref,  7>  we  hypothesize  that  the  film  thickness 
can  be  expressed  in  the  form 

h  =  h*  exp£^(x)/  (A7) 

where  ^.(x)  and  its  derivatives  are  of  0(1). 

■dp  ^ 

Nov;  let  p  and  v  be  the  fluid  density  and  viscosity, 
respectively,  corresponding  to  the  state  (p*,T*).  Then  de¬ 
fine  the  new  dependent  variables: 


U.  — 


u  L 

' 


(A8) 


Other  barred  quantities  are  to  be  made  by  dividing  the  local 
dimensional  quantity  by  its  value  at  the  reference  state 
(p*,T*).  Thus:  k  =  kA(p*,T*),  etc. 

The  barred  variables  are  now  substituted  into  eqs.  1-4. 
When  these  substitutions  are  made,  the  characteristic  small 
parameter  of  fluid-dynamic  lubrication  appears;  i.  e., 

S  =  h*/L  (A9) 


The  resulting  equations  are: 


2*7- 


-£(¥4 

+  (A10) 


(All) 


40 


THE  FRANKLIN  INSTITUTE 


3-fl  (?*)+£%?*>  =° 


Laboratories  for  Research  and  Development 

I-A20h9"10 

(A12) 


(A13) 

The  boundary  conditions  to  be  satisfied  are: 
u  =  UL/P*,  7  =  0,  f  =  1;  7=0 

u  =  0  ,  v  =  0,  f  ~  1;  y=l  (Al4) 

P  =  1;  5c  =  0  or  x  =  1 

Expansion  in  b  : 

As  in  ref.  7,  let  us  now  restrict  consideration  to  those 
solutions  of  the  system  10-14  which  can  be  expanded  in  US  ", 
Thus,  we  take: 

d  Ug  etc. 


U  =  U0  + 


v  =  V0  +  Sv1  + 


(A15) 


p  =  p0  +  +  . 

I  =  T0  +  5lj_  +  . 

Here  the  capitalized  functions  are  to  be  independent  of  " <5  n. 
The  zeroth-order  functions  are  made  to  assume  all  of  the  non¬ 
zero  boundary  conditions  for  their  respective  barred  quanti¬ 
ties  so  that  all  higher  order  functions  must  vanish  on  the 
boundaries. 
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Before  the  substitution  of  functions  .15  into  eqs.  10-14 
is  made,  one  should  note  that  all  x  derivatives  at  constant 
" y "  can  be  converted  to  3T  derivatives  at  constant  with¬ 

out  involvement  of  "  S  " .  Thus  if  ^  is  any  function  of 
position,  v/e  have: 


All  functions  of  the  thermodynamic  state  can  be  assumed 
to  deviate  from  the  state  (PQ,  T  )  by  0(«£).  Thus: 


-/>  (z,  v  +  ^0Jr7/'  +  ^(W)p  ^  (A2o) 

When,  now,  the  series  15  and  20  are  substituted  into  eqs.  10-14, 
and  the  coefficients  of  the  various  powers  of  " £  "  are  equated 


to  zero,  we  get: 


(A21) 

(A22) 


(A23) 
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Similar,  though  more  complicate!,  differential  equations 
are  obtainable  for  the  higher  order  functions,  The  equations 
for  laminar  lubrication  are  obtained  bv  approximating  the 
barred  quantities  in  eq.  15  by  their  zeroth  order  functions* 
When  reversion  is  made  to  the  original  physical  variables, 

t\  h 

and  allowance  is  made  for  the  additional  space  variable,  z, 
the  resulting  equations  are  eqs.  1-5  of  the  main  text. 


< 


43 


THE  FRANKLIN  INSTITUTE  •  Laboratories  for  Research  and  Development 


I- A2049- 10 


DISTRIBUTION  LIST 


No.  of  Copies 


No.  of  Copies 


10 

X 

1 

1 

1 

1 


Chief  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  DC. 
Attn:  Code  429 
Code  461 
Code  423 
Code  438 
Code  463 
Code  466 


1 


1 


1  Commanding  Officer 

Office  of  Naval  Research 

Branch  Office  1 

346  Broadway 

New  York  13,  New  York 

1  Commanding  Officer 

Office  of  Naval  Research 

Branch  Office  1 

1030  East  Green  Street 

Pasadena  1,  California 

1  Commanding  Officer 

Office  Naval  Research 

Branch  Office  3 

Tenth  Floor 

The  John  Crerar  library  Building 
86  East  Randolph  Street 
Chicago,  Illinois 


Commander 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base. 

Ohio 

Attn:  WCLJN  -  Prof.  Wm,  J.  Moreland 
Commander 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base. 

Ohio 

Attn:  WCLP'.-l  -  Mr.  G.  L.  Wander 
Commander 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base, 

Ohio 

Attn:  WCLPF-2  -  Mr.  G. A.  Beane 
Commander 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base, 

Ohio 

Attn:  WCLSM-  Mr.  P.C.  Hanlon 
Commander 

Air  Force  Office  of 
Scientific  Research 
Washington  25,  D.C. 

Attn:  Mr.  Joseph  E.  Long  (SRI1A) 


10  Commanding  Officer  3 

Office  of  Naval  Research 
Branch  Office 

Navy  #100  Fleet  Post  Office 
New  York,  New  York 

8 

1  Office  of  Naval  Research 

Resident  Representative 
University  of  Pennsylvania 
Room  213,  Hare  Building 
Philadelphia  4,  Pennsylvania  2 

1  Director 

Naval  Training  Device  Center 
Sands  Point,  Port  Washington 
Long  Island,  New  York  4 

Commanding  General 
US  Army  Engineer  Research  and 
Development  Laboratories 
Fort.  Belviir,  Virginia 

2  Attn:  Mr,  W.M.  Crim,  Nuclear 

Power  Field  Office 

1  Technical  Intelligence  1 

Branch  6 


Chief,  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.C. 

Attn:  Mr.  Norman  (AE  6333) 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.C. 

Attn:  Mr.  James  C.  Reid  (Code  644B) 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.C. 

Attn:  Mr.  0. J.  Stanton  (ReS-1) 

Chief,  Division  of  Engineering 
Maritime  Administration 
GAO  Building 
Washington  25,  D.C. 

Di  rector 

Nava]  Research  Laboratory 
Washington  25,  D.C. 

Attn:  Code  5230 
Code  2000 


Mr.  Norman  L.  Klein 
Research  and  Development  Division 
Office  of  Chief  of  Ordnance 
Department  of  the  Army 
Washington  25,  D.C. 


Director 

US  Naval  Engineering  Experiment  Station 
Annapolis,  Maryland 
Attn:  Supervisor,  Bearings  Project, 
Code  790 


1  Commander 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base, 
Ohio 

Attn:  WCLJN  Mr.  J.C.  Noble 


Army  Reactors  Branch 
Reactor  Development  Division 
US  Atomic  Energy  Commission 
Washington  25,  D.C. 

Attn:  Mr.  Clarence  Miller 
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1  Chief,  Aircraft  Reactor  Branch 

Division  of  Reactor  Development 
US  Atomic  Energy  Commission 
Washington  25,  D.C. 

Attn:  Maj .  Gen.  D.J.  Keirn 

3  Technical  Library 

US  Atomic  Energy  Commission 
Washington  25,  D.C. 

Attn:  John  L.  Cook 

1  Jet  Propulsion  Laboratory 

California  Institute  of  Technology 
4800  Oak  Grove  Avenue 
Pasadena,  California 
Attn:  Dr.  John  H.  Laub 

1  Chief,  Technical  Information 

+■  repro  Service  Extension 
P.0.  Box  62 
Oak  Ridge,  Tennessee 
Attn:  Mr.  Melvin  S.  Day 

1  Head,  Experimental  Engineering,  R.P.D. 
Oak  Ridge  National  Laboratory 
Building  9201-3 

Oak  Ridge,  Tennessee 
Attn:  Mr.  H.W.  Savage 

2  Director 

National  Aeronautics  and 
Space  Administration 
1512  H.  Street,  N.W. 

Washington  25,  D.C. 

Attn:  Mr,  Harold  Hessing 

Chief  of  Staff,  US  Air  Force 
The  Pentagon 
Washington  25,  D.C. 

1  Attn:  DCS/D  AFDRD/AD-2 

1  DCS/D  AFOP-OO 

1  Chief  of  Research  and  Development 

Office,  Chief  of  Staff 
Department  of  the  Army 
Pentagon  Building 
Washington  25,  D.C. 

1  Commanding  Officer 

Diamond  Ordnance  Fuze  Laboratories 
Washington  25,  D.C. 

Attn:  Technical  Reference  Section 
(CRDTL  06,33) 

1  Office,  Chief  of  Engineers 

Engineer  Research  and 
Development  Division 
Gravelly  Point 
Washington  25,  D.C. 

1  Commanding  Officer 

Office  of  Ordnance  Research 
Box  CM,  Duke  Station 
Durham,  North  Carolina 


1  Commanding  Officer 

Detroit.  Arsenal 
CRD  MX-ECPD 
Centerline,  Michigan 

1  Commander, 

Army  Rocket  and  Guided  Missile  Agency 

Redstone  Arsenal 

Alabama 

Attn;  Technical  Library,  QRDXR-OTL 

1  Commanding  Officer 

US  Naval  Boiler  and  Turbine 

Laboratory 

Naval  Base 

Philadelphia,  Pennsylvania 

1  Commanding  Officer 

US  Naval  Avionics  Facility 
Indianapolis  18,  Indiana 
Attn:  Mr.  J.  G.  Weir 

1  Superintendent 

US  Naval  Postgraduate  School 
Monterey,  California 

1  Commander 

Air  Research  and  Development  Command 
Andrews  Air  Force  Base 
Washington  25,  D.C. 

Attn:  Lt.  Col.  C.D.  Reifstech, 

RETAP 

1  Committee  on  Equipment  and  Supplies 

Office  of  the  Assistant  Secretary 
of  Defense  (R&D) 

The  Pentagon 
Washington  25,  D.C. 

Attn:  Mr.  A.F.  Bird 

1  Technical  Library 

Office  of  the  Assistant 
Secretary  of  Defense  (R&D) 

The  Pentagon 
Washington  25,  D.C. 

Attn:  Mr.  William  H.  Plant 

10  Armed  Services  Technical 

Information  Agency 
Arlington  Hall  Station 
Arlington  12,  Virginia 
Attn:  TIPCR 

1  Office  of  Technical  Services 

Department  of  Commerce 
Washington  25,  D.C, 

1  US  Department  of  Commerce 

National  Bureau  of  Standards 
Boulder  Laboratories 
Boulder,  Colorado 

1  A.C.  Spark  Plug  Division 

General  Motors  Corporation 
Milwaukee  1,  Wisconsin 
Attn:  Mr.  Allen  Knudsen 
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1  Aerojet-General  Corporation 

Nucleonics 

San  Ramon,  California 

1  AiResearch  Manufacturing  Division 

The  Garrett  Corporation 
9851  S  Sepulveda  Boulevard 
Los  Angeles,  California 
Attn:  Dr.  W.T.  Von  Der  Nuell 

1  Allis  Chalmers  Manufacturing  Co. 

Buda  Division 
1126  South  70th  Street 
Milwaukee  1,  Wisconsin 
Attn:  Mr.  H.K.  Ihrih 

Vice  President  for  Research 

1  American  Society  of  Lubrication 

Engineers 

5  North  Wabash  Avenue 
Chicago  2,  Illinois 

1  Chairman,  Research  Committee  on 

Lubrication 

American  Society  of  Mechanical  Fngi 
29  West  39th  Street 
New  York  18,  New  York 

1  Analogue  Controls,  Incorporated 

200  Frank  Road 
Hicksville,  L  I  ,  New  York 
Attn:  Mr.  J.  L.  Cherubin, 

Chief  Engineer 

1  Applied  Physics  Laboratory 

Johns  Hopkins  University 
Silver  Spring,  Maryland 
Attn:  George  L.  Seielstad, 
Supervisor,  Technical 
Reports  Group 

1  Autonetics 

A  Division  of  North  American 
Aviation,  Inc. 

9150  East  Imperial  Highway 
Downey,  California 

1  The  Barden  Corporation 

Danbury,  Connecticut 
Attn:  Mr.  B.L.  Mims,  Vice  Pres. 

1  Battelle  Memorial  Institute 

505  King  Avenue 
Columbus  1,  Ohio 
Attn:  Dr.  Russell  Dayton 

1  Bell  Aircraft  Corporation 

P.O.  Box  1 
Buffalo  5,  New  York 
Attn:  Mr.  Roy  J.  Sand strom 
Vice  President  (Engr) 


No.  of  Copies 

1  'Jtica  Division 

Bendix  Aviation  Corporation 

Utica,  New  York 

Attn:  Mr.  Russell  T  DeM-:"h, 

Senior  Engineer 

1  Boeing  Airplane  Company 

Research  Group 
Box  3707 

Seattle  25,  Washington 
Attn:  Mr.  D.J.  Melchior 

Mechanical  Equipment 

1  Bryant  Chucking  Grinder  Company 

60  Clinton  Avenue 
Springfield,  Vermont 
Attn:  Mr.  Hugh  Taft 

1  Cadillac  Gage  Company 

P.O.  Box  3806 
Detroit  5,  Michigan 
Attn:  Mr.  J.  Taylor, 

Project  Engineer 

1  Carrier  Corporation 

Division  of  Research 
Research  Center 
Syracuse,  New  York 
Attn:  Dr.  Dewey  J.  Sandell 

1  Chance  Vought  Corporation 

P.O.  Box  5907 
Dallas,  Texas 
Attn:  Mr.  R.C.  Blaylock 

Vice  President  (Engr) 

1  Chrysler  Corporation 

Defense  Operations 
P.O.  Box  757 
Detroit  31,  Michigan 
Attn:  Mr.  C.W.  Snider 

1  The  Cleveland  Graphite  Bronze  Co. 

17000  St.  Clair  Avenue 
Cleveland  10,  Ohio 
Attn:  Mr.  R.  H.  Josephson 

1  Convair 

A  Division  of  General  Dynamics  Corp. 

Fort  Worth,  Texas 

Attn:  Mr.  L.E.  McTaggart 

Senior  Design  Engineer 
Engineering  Annex  #1 

1  Wright  Aeronautical  Division 

Curtiss-Wright  Corporation 
Department  8311 
Wood -Ridge,  New  Jersey 
Attn:  Mr.  W.  J.  Derner 

Chief  Project  Engineer 
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1  Daystrom  Pacific 

9320  Lincoln  Boulevard 
Los  Angeles  45.  California 
Attn;  Mr.  Don  Ensley 

Special  Projects  Engineer 

1  Detroit  Controls  Division 

American  Radiator  &  Standard 
Sanitary  Corp. 

100  Morse  Street 
Norwood.  Massachusetts 
Attn;  Mr.  E.L.  Swainson, 

Chief  Engineer 

1  E. I.  duPont  de  Nemours  &  Co. , Inc. 

Wilmington  98,  Delaware 
Attn:  Dr.  McNeilly,  Mechanical 
Development  Laboratory 

l  Dynamic  Controls  Corporation 

14  Runde  Lane 
Bloomfield,  Connecticut 
Attn:  Mr.  Thomas  P.  Farkas 

1  Electronics  Systems  Laboratory 

Federal  Telecommunications  Labs. 
Nut  ley.  New  Jersey 
Attn:  Mr.  John  Metzger 

1  Ford  Instrument  Company 

31-10  Thomson  Avenue 
Long  Island  City  1,  New  York 
Attn:  Mr.  Jarvis 

1  Ford  Motor  Company 

Dearborn.  Michigan 
Attn:  Professor  Adolph  Egli 

1  General  Atomics  Division 
General  Dynamics  Corporation 
P.0.  Box  608 

San  Diego  12.  California 
Attn:  Mr.  F.W.  Simpson 

2  General  Engineering  Laboratory 
General  Electric  Company 

One  River  Road 
Schenectady  5,  New  York 
Attn:  Dr.  Beno  Sternlicht 

1  General  Motors  Corporation 

Research  Laboratories  Division 
Technical  Center 
Detroit  2,  Michigan 
Attn:  Mr.  Robert  Davies. 

.Assistant  Supervisor 
of  Bearings  Section 

1  Grumman  Aircraft  Engineering  Corp 

Bethpage,  Long  Island,  New  York 
Attn:  Mr.  John  Karanik,  Chief 

Power  Plant  &  Equipment 
Installation  Department 


No.  of  Copies 

1  Institute  of  Aeronautical  Sciences 

2  East  64th  Street 
New  York,  New  York 

1  International  Business  Machine  Corp. 

Research  Laboratory 
San  Jose,  California 
Attn:  Dr.  W.  A.  Gross 

1  International  Tel.  &  Tel  Labs. 

15151  Bledsoe  Street 
San  Fernando,  California 
Attn:  Mr.  G.  B.  Speen 

1  Jack  &  Heintz,  Incorporated 

207  East  Laboratory 
Public  Roads  Building 
Washington  National  Airport 
Washington  1,  D.C. 

1  Kearfott  Company 

1378  Main  Street 
Clifton,  New  Jersey 
Attn:  Mr.  Walter  Carow 

l  James  J.  Barker 

Associate  Professor  of  Nuclear  Engrg 
Department  of  Chemical  Engineering 
New  York  University 
New  York  53.  New  York 

1  Koppers  Company,  Incorporated 

Metal  Products  Division  > 

P.0.  Box  626 
Baltimore  3,  Maryland 
Attn:  Mr.  T.  C.  Kuchler 

1  Lear,  Incorporated 

Grand  Rapids  Division 
110  Ionia  Avenue,  N.W. 

Grand  Rapids  2,  Michigan 

1  Lear,  Incorporated 

3171  South  Bundy  Drive 
Santa  Monica,  California 
Attn:  Richard  M.  Mock,  President 

1  Litton  Industries 

336  North  Foothill  Road 
Beverly  Rills,  California 
Attn:  Mr.  D.  Moors,  Department  21 

1  Lockheed  Aircraft  Corp. 

Missiles  and  Space  Division 
P.O.  Box  504 
Sunnyvale,  California 
Attn:  Dr.  M.  A.  Steinberg,  53-30 
Bldg.  202.  Plant  2,  PA 

1  Lycoming  Division 

Avco  Manufacturing  Corporation 
Main  Street 
Stratford,  Connecticut 
Attn:  Mr.  S.  B.  Withington 

President  &  General  Manager 
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McDonnell  Aircraft  Corporation  1 

Lambert  St.  -  St.  Louis 
Municipal  Airport 
Box  516 

St.  Louis  3,  Missouri 
Attn:  Mr.  Kendall  Perkins 

Vice  President  (Engr.)  1 

Marquardt  Aircraft  Company 
P.0-  Box  2013  -  South  Annex 
Van  Nuys.  California 

Attn:  Mr.  A.J.  Kreiner,  Director  1 

Controls  &  Accessories  Div. 

Massachusetts  Institute  of  Technology 
Department  of  Aeronautical  Engineering 
Instrumentation  Laboratory 
Cambridge  39,  Massachusetts  1 

Attn:  Dr.  L  Mel).  Schetky 

Dynamic  Analysis  and  Control  Lab. 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  1 

Attn:  Dr.  R.  W.  Mann 

Aeronautical  Engineering  Laboratories 
University  of  Michigan 
Ann  Arbor,  Michigan 

Attn:  Professor  R.  B.  Morrison  1 

Minneapolis-Ilonywel  1  Regulator  Co. 

Gyro  Section 

2600  Ridgeway  Road  1 

Minneapolis,  Minnesota 
Attn:  Mr.  Jack  W.  Lower 

Chief,  Aeronautical  Division 


Ryan  Aeronautical  Company 
Lindberg  Field 
San  Diego  12,  California 
Attn:  Mr.  Frank  W.  Fink 

Vice  President  &  Chief  Engr. 

Sanderson  S  Porter 
72  Wall  Street 
New  York  5,  Nev.-  York 
Attn-  Mr.  S.  T.  Robin  son 


Smith  Bearing  Company 

Subsidiary  of  the  Ex-Cello-0  Corn. 

P.0  Box  1119 

Trenton,  New  Jersey 

Attn:  Mr.  L.  A.  Mass 


Solar  Aircraft  Companv 
2200  Pacific  Highway 
San  Diego  12,  California 
Attn:  Paul  A.  Titt,  Chief  Engr. 

Sperry  Gyroscope  Company 
Division  of  Sperry  Rand  Company 
Great  Neck,  Long  Island,  N.Y. 
Attn:  Mr.  Wilmer  L*  Barrow 

Vice  President  for  R&D 


Stevens  Institute  of 
Hoboken,  New  Jersey 
Attn:  Professor  P.F. 


Technology 

Martinuzzi 


Stratos  Division 
Fairchild  Airplane  &  Engine  Co. 
Hay  Shore,  Long  Island,  New  York 
Attn:  Mr.  George  Mackowski 


New  Departure  1 

Division  of  General  Motors  Corporation 
Bristol,  Connecticut 
Attn:  R.N.S.  Scheidt,  Research 

Test  Laboratories  Supervisor  I 

Northrop  Aircraft,  Incorporated 
Northrop  Field 

Hawthorne,  California  1 

Attn:  Mr.  Robert  R.  Miller 

Vice  President  &  General  Manager 

The  Technological  Institute 

Northwestern  University  ^ 

Evanston,  Illinois 

Attn:  Professor  A.  Charnes 


Sundstrand  Turbo 
2480  West  70th  Avenue 
Denver  21,  Colorado 


Sundstrand  Turbo 
10445  Glenoaks  Boulevard 
Pacoima,  California 


ihompson  Ramo  Wooldridge  Corporation 
P.O.  Box  90534  Airport  Station 
Los  Angeles  45,  California 
Attn:  Mr.  Robert  I.  King  F1829 


Tribo-Netics  Laboratories 

Vermilion,  Ohio 

Attn:  Mr.  E.  Fred  Macks 


Owens-Corning  Fiberglass  Corporation 
Pioneering  Laboratory 
P.O.  Box  415 
Granville,  Ohio 

Attn:  Mr.  R.L.  Boithwick,  Project  Manager 


sxi,  tenter 

3901  N.E,  12th  Avenue 
Pompano  Beach,  Florida 
Attn:  Dr  W.  C.  Knopf 


1  Badio  Corporation  of  America 

Camden,  New  Jersey 
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1  Westinghouse  Electric  Corporation 

Research  Laboratories 
East  Pittsburg,  Pennsylvania 
Attn:  Mr.  John  Boyd. 


1  Westinghouse  Electric  Corporation 

P.O.  Box  868 

Pittsburg  30,  Pennsylvania 
Attn:  Mr.  J.  K.  Hardnette 

Vice  President  and  General  Man. 


1  Worthington  Corporation 

Harrison,  New  Jersey 
Attn:  Mr,  H.  Walter 

Director  of  Research 


1  Universal  Match  Company 

Avionics  Department  Technical  Library 

4407  Cook  Avenue 

St.  Louis  13,  Missouri 


No.  of  Copies 

1  Maico  Electronics  Inc. 

21  North  Third  Street 
Minneapolis  1,  Minnesota 
Attn:  Mr.  G  E.  Adams, 

Automatic  Controls  Dept. 

1  Collins  Construction  Company 

Post  Office  Box  86 
Port  Lavaca,  Texas 
Attn:  Mr.  James  11.  Mamarchev 
Consulting  Engineer 

1  Space  Technology  Laboratories,  Inc. 

P.O.  Box  95001 
Los  Angeles  45,  California 
Attn:  Mr.  Kenton  Andrews 

Building  "F1',  Boom  125 

1  Beemer  Engineering  Company 

401  North  Broad  Street 
Philadelphia  8,  Penna. 

New  Hampshire  Ball  Bearings  Inc 
Peterborough,  New  Hampshire 
Attn:  Mr.  Henry  F.  Villaume,  Engr. 
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